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Abstract 

We study a class of nonperturbative corrections to single-inclusive photon cross sections 
at measured transverse momentum pr, in the large-XT limit. We develop an extension of 
the joint (threshold and transverse momentum) resummation formalism, appropriate for 
large xt, in which there are no kinematic singularities associated with recoil, and for which 
matching to fixed order and to threshold resummation at next-to-leading logarithm (NLL) 
is straightforward. Beyond NLL, we find contributions that can be attributed to recoil from 
initial state radiation. Associated power corrections occur as inverse powers of p\ and are 
identified from the infrared structure of integrals over the running coupling. They have 
significant energy dependence and decrease from typical fixed-target to collider energies. 
Energy conservation, which is incorporated into joint resummation, moderates the effects 
of perturbative recoil and power corrections for large xt- 
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1 Introduction 



To study the interplay of perturbative and nonperturbative dynamics in processes involving 
hadronic states it is natural to begin with observables whose perturbative analysis is well un- 
derstood. For certain observables, perturbation theory not only provides predictions at leading 
power in a large momentum scale, but also characterizes power corrections in that scale. This can 
come about, for example, through nonconvergent perturbative expansions that exhibit sensitivity 
to the strong- coupling and/or vacuum structure of the theory [1]. 

Relying on perturbative resummations, this approach has had phenomenological successes in 
the description of a variety of inclusive and semi-inclusive cross sections. These include average 
and differential event shapes, primarily but not exclusively in e + e~ annihilation [2-9], and elec- 
troweak annihilation cross sections at measured transverse momentum [10-13]. The value of the 
event shape or lepton pair transverse momentum provides a second scale in the cross section, and 
varying this scale changes the relative importance of perturbative and nonperturbative dynamics. 
Thus, the transition between perturbative and nonperturbative QCD is in principle available for 
study in these observables. 

In this paper, we adopt this general philosophy and employ the joint resummation [14, 15] 
of threshold [16] and transverse momentum [12] enhancements to study power corrections in the 
hard-scattering scale px for single-particle inclusive (1PI) cross sections in the large x\ = kp\jS 
region. Using direct photon production as an example, we will show that these corrections exhibit 
significant xt dependence, which moderates both perturbative and nonperturbative recoil at large 
xt compared to estimates based on transverse momentum resummation alone. These conclusions 
are made possible by a simplification of the joint resummation formalism that is specific to the 
Xt — > 1 limit. 

Direct photon production was originally envisioned as a relatively straightforward process 
with which to test fixed-order perturbative calculations and to determine the gluon distribution 
[17, 18, 19]. The extensive data on direct photon production [20-24], however, has turned out 
to be more complex than was perhaps expected. Presumably for this reason it has inspired 
varied theoretical and phenomenological studies [25-37]. Nevertheless, for this benchmark process 
important questions remain unresolved. In particular, it has been argued that fixed target data 
for direct photon production in the lower pt range (roughly below 5 GeV) are difficult to reconcile 
with collinear-factorized NLO cross sections [29, 31]. Additionally, threshold resummation [33, 
34, 35] appears to explain the data only for larger p?. 

This difficulty has motivated the use of /^-dependent, or unintegrated, parton distributions 
combined with recoiling partonic 2 — * 2 subprocesses [30] . Information on the partonic transverse 
momenta in such distributions may come from resummed perturbation theory [12, 13, 30], and/or 
from comparisons to data [11, 29, 38, 39], including Drell-Yan, photon and hadron pair cross 
sections. Probably the simplest approach is to assume a Gaussian dependence exp[— k T /(k T )] 
[18]. As we review below, perturbative resummations predict logarithmic pt dependence for the 
parameter (k T ). They also imply that (k T ) depends on the parton flavor. 

The use of unintegrated distributions requires an extension of collinear factorization [40] . In 
particular, a technical challenge in the case of light particle production is the potential for an 
artificial infrared singularity when the total transverse momentum of the initial state partons 
is comparable to the observed p? [18]. One way to avoid this singularity is to impose strong 
ordering in transverse momenta, as in [30] , a procedure which requires definition beyond leading 
logarithm. A related approach is described detail in Ref . [32] , based on a specific implementation 
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of /cr-resummation. As presented in [32], however, fits in this formalism favor fixed (kj), with 
no indication of the pr-dependence implied by kx resummation. Other studies, however, seem to 
imply that (k^) is S'-dependent [29]. In summary, it remains unclear how much of what we inter- 
pret as recoil, or parton transverse momentum, is perturbative and how much nonperturbative. 
Here we come back to this question in the context of a generalized resummation formalism. 

Resummed perturbation theory for 1PI cross sections was extended in [14, 15, 37] using joint 
resummation. Joint resummation systematically combines singular behavior at zero transverse 
momentum for initial-state partons with that at partonic threshold, where the initial state par- 
tonic invariant mass s = x a XbS is just large enough to produce the observed final state. This 
method was applied to Z and Higgs production in [41, 42], where no kinematic singularities arise, 
because the produced electroweak state is massive and the transverse momentum of the lepton 
pair is directly observed. In [41], some implications for the specific forms of power corrections 
were also pointed out. Although the joint formalism was applied to high-p-r photon production in 
[37], its application was hampered by the same infrared singularity mentioned above, associated 
with the production of a massless particle. As noted in [37], the complexity of a simultaneous 
resummation in transverse momentum and energy above threshold appears to make impractical 
a matching of the sort developed for transverse momentum resummation alone in [32]. 

In this paper, we extend this work, and revisit logarithmic and power corrections to the direct 
photon cross section in the joint resummation formalism. Compared to previous work, however, 
we use the kinematics of the large-x^ limit to reformulate joint resummation, taking into account 
recoil effects in the partonic subprocess while avoiding a kinematic singularity. The resulting 
resummed cross section reduces to threshold resummation at next-to-leading logarithm (NLL) 
and can be matched to finite-order and threshold resummed cross sections in a straightforward 
fashion. Beginning at NNLL, the cross section also includes a contribution that can be identified 
as the finite residue left from the cancellation of the transverse momentum singularities of real 
and virtual gluons radiated in the initial state. Enhancements to the cross section associated 
with final state interactions are treated only to leading power in this paper, and appear in the 
same manner as in threshold resummation. In another paper we will argue that the results found 
here are not changed qualitatively by these effects. 

The parameters that control power corrections associated with joint resummation at par- 
tonic threshold are found to be related to parameters familiar from the transverse momentum 
distributions in electroweak annihilation. The power corrections also inherit significant energy 
dependence. For large Xt, both perturbative recoil and nonperturbative power corrections to 
the predictions of threshold resummation are suppressed by the phase space restrictions built 
into joint resummation. This effect is important, however, only for xt near one, or equivalently 
for large values of its conjugate Mellin moment variable N. For smaller xt or N of order unity 
an analysis based on kr resummation alone may be appropriate, but should be matched to the 
results of joint resummation in the large xt region. 

We begin Sec. 2 with a brief summary of the joint resummation formula as developed for direct 
photon cross sections, and exhibit the kinematic singularity. In the next subsection, the cross 
section is expressed as a double inverse transform. This is followed by a simple reformulation 
that eliminates the kinematic singularity and reduces the jointly resummed cross section to a 
single transform that extends threshold resummation for the direct photon cross section. The 
resulting Sudakov exponents of joint resummation are analyzed in Sec. 3, where we identify the 
form of the recoil and power corrections to the direct photon cross section that are implied by 
joint resummation. We explore the phenomenology of these corrections in Sec. 4, exhibit the 
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suppression of power corrections for large S at fixed pt, and briefly discuss possible subdominant 
corrections not directly associated with partonic threshold. We conclude with a summary, and 
a brief discussion of possible implications for an eventual global treatment of single-photon and 
single-hadron cross sections. 

2 Self-consistent Recoil in Joint Resummation 
2.1 Partonic recoil in direct photon production 

Joint resummation [14, 37] is an extension of threshold [16] and transverse momentum resum- 
mations [12, 13] that unifies these two formalisms. So far, at the phenomenological level it has 
been applied primarily to the single electroweak boson (mass Q) production cross sections at low 
transverse momentum, Qt <Q [41, 42]. In this case, threshold resummation is associated with 
corrections of the form [a s n /(l — z)] ln 2ra_1 (l — z), with z = Q 2 /s, where y/l is the invariant mass 
of the partonic pair that annihilates into the observed boson. Such corrections are "implicit" 
in the sense that they contribute to the hadronic cross section only after convolution with the 
parton distribution functions, and hence give nonlogarithmic, although potentially significant, 
contributions to the cross section. Singular corrections in Qt, on the other hand, are explicit 
in the cross sections themselves, appearing as terms like a s n /Q 2 T \n 2n ~ l (Q T /Q) directly for the 
measured spectrum. 

For single-particle inclusive cross sections such as direct photon production at measured 
p T 3> Aq CD the situation is slightly different. To leading order in the hard scattering, incoming 
partons produce a photon-parton system, which subsequently evolves into a photon-jet pair. At 
higher orders in a s , the pair recoils against unobserved soft gluon radiation with total transverse 
momentum Qt, in much the same way as for a single electroweak boson. When only the photon 
is observed, Qt is integrated and singularities at Qt/pt = cancel, analogously to singularities 
at 1 — z — in threshold resummation. Thus in the direct photon cross section at measured 
Pt, both transverse momentum and threshold singularities are implicit rather than explicit. 
Nevertheless, small-Q^ gluon radiation can play a significant role in the cross section. Applied 
to direct photon production, joint resummation attempts to estimate the effects of these soft 
emissions systematically [14, 37]. 

For an observed photon of momentum pt, the photon transverse momentum in the pair 
center-of-mass is 

Pt = Pt~Qt/2. (1) 

In the limit that Qt/pt *C 1 the cross section is a convolution [37] of the resummed distribution in 
Qt with a hard-scattering function evaluated at photon momentum p' T . Non-zero pair momentum 
Qt, if in the direction of the observed photon, decreases the scale of the hard scattering, and 
can thus enhance the cross section. As emphasized in Ref. [37], however, when Qt grows to the 
order of pt, this approximation generates kinematic singularities. Their effect is non-negligible 
because the fall-off in soft gluon transverse momenta has a power-like perturbative tail. In [37], 
we dealt with the kinematic singularity in a rather crude way by cutting off the resummed Qt 
spectrum at a convenient scale Qt = fi'- 

J res r- rlrr vcs 
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where Q is the invariant mass of the photon-parton pair. At threshold, the latter is fixed by 

Q = 2p T . (3) 

The scale p, in Eq. (2) may be regarded as a matching scale. Ideally, at Qt ~ p> < Pt, one 
would replace the resummed cross section p%da r X s B -> 1 x I ^Q 2 d 2 Qr dpr by the fixed-order (NLO) 
one, which does not have the kinematic singularity. In practice, this becomes a very complicated 
procedure, and it is more convenient to derive a jointly resummed cross section that does not 
require a cutoff. We will show below that this may be achieved by applying an additional, self- 
consistent approximation that is exact at partonic threshold. To do so, we must recall the explicit 
form of the cross section derived in Refs. [14, 37]. 



2.2 The double inverse transform 

Integrated over rapidities, the jointly resummed direct photon cross section is written in terms 
of Mellin moments of the MS parton distributions, 4> a /n{N, fi) = L dxx N ~ Y </> a /#(x,yu), as 



3 da { ™] x v \ [dN~ ~ f 1 2 , 2 . N \M ab (x 2 T )\ 2 

x C ab ^ c (a s (p),x 2 T ) J^®(P-Qt) ^ P ab (iV,Q r ,^,^ , (4) 



where /i is the factorization and renormalization scale, and the |M a b| 2 are squared amplitudes for 
the partonic processes ab —> 7c. The variable defined by 

cosh T) 

where fj is the rapidity of the direct photon in the center of mass of the hard scattering. At 
partonic threshold, or equivalently large values of the moment variable N, fj is forced to unity. 
For this reason, we will approximate 

2p T /x T ~ 2p T = Q (6) 

in the functions P a b in Eq. (4), where dependence on pt is logarithmic. The contour C in Eq. (4) 
and the b integral in (7) below define the inverse transforms from N, b space to z and Qt- These 
contour integrals were described in detail in Refs. [37, 41]. 

The functions P a b in Eq. (4) were derived in Ref. [37] and provide Qt dependence at fixed N. 
Each P a b is itself the Fourier transform of the exponentiated logarithmic dependence on N and 
b, 

P ab (N, Q T , Q,ijl) = J d 2 be~ ib - QT exp [E ab ^ c (N, b, Q, /x)] , (7) 

where the E a b^ lc are "Sudakov" exponents that we will specify explicitly below. They can be 
split into initial and final state contributions, where, as shown in [37], all 6-dependence comes 
from the initial state, 

E ab ^ c (N, b, Q, n) = EH(N, b, Q, 11) + E F ab s c (N, Q, fx) . (8) 
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The iV-independent coefficients C ^ c contain the effects of hard virtual corrections and are 
perturbative series of the form C ah ~^ yc = 1 + ^.(7 ctfe - > 7c(i) _)_.,,. To next-to-leading logarithmic 
accuracy one needs the first order terms which may be found in [34], and are given in Appendix 
A below. In this paper, we concentrate on the initial state exponent, which contains all leading 
logarithmic effects and all h dependence. 

Threshold resummation is recovered from Eq. (4) by setting b to zero in the exponents E a ^ lc . 
In this case, the b integral produces a delta function that sets p' T = pt- Then the exponent 
E a b^ c reverts to its threshold resummed form, and S/4p' T 2 — > \jx\. Recoil enhances the cross 
section Eq. (4) because even for Q T <C p T , the ratio S/Ap' T 2 can be larger than S/A.p T . For 
large enough Q T , the factor S/4p' T 2 can diverge, and a cutoff is required, as discussed above. 
This momentum configuration, however, requires, Qt ~ %Pt, and hence is far outside the region 
where resummation is applicable. This problem is not due to our approximation in the region of 
interest, but to our extrapolation beyond that region. When Q T approaches px in magnitude, 
the factorization between gluon emission and hard-scattering fails. For large N, however, the 
"profile" functions P a b vanish once Qt > Q/N <C pr- The P a bS vanish for moderate Qt because 
the exponents E a b_ lc develop large (negative) logarithms once bQ/N > 1 (see below, Eq. (17)). 
This ensures that for Qt > Q/N the exponential exp[— ib ■ Qt] oscillates on a smaller scale than 
the size of the region where exp[i? a (,_> 7C ] is nonvanishing. Put another way, because the widths 
of the profile functions in b space are of order N/Q, their Fourier transforms to Qt space have 
widths of order Q/N. Numerical examples for the Qt integrand in Eq. (2) were given in Ref. 
[37], which show the fall-off of the profile function for increasing Qt, followed by the kinematic 
singularity as Qt increases to the order of pt- We conclude that for large N the true enhancement 
due to recoil is insensitive to modifications of the integrand above Qt ~ Pt/N <C pt- Since large 
N corresponds to xt 1, we expect a suppression of recoil effects in this limit. In addition, 
N is conjugate to /co/pt, where ko is the energy of initial state radiation [37, 14]. The relation 
Qt < Pt/N is thus equivalent to the restriction that the total transverse momentum of initial 
state radiation is less than its energy. We will use this observation shortly. 

2.3 Elimination of the kinematic singularity 

Given that all-order recoil effects enhance the jointly resummed cross section from values of 
Qt such that Qt < Pt/N, it is only in this region that we are required to maintain accurate 
expressions for leading Qt behavior (that is, 1/Q T times logarithms). In fact, to construct the 
jointly resummed expression in Eq. (4), we have neglected corrections that are nonsingular at 
Qt = and 1 — z — 0. This means that we do not in general have control over corrections 
suppressed by powers of Qt/pt, and also that we are free to change the resummed expression 
at this level of accuracy. Such a change will only affect the result from the region of Qt beyond 
the range that gives enhancement. These modifications will not produce logarithms, and we can 
adjust for them by matching to the cross section at fixed order. 

In summary, we are free to choose an extrapolation that does not produce spurious singu- 
larities at large Qt and which does not change the singularity structure at Qt = 0. This may 
be done in such a way that the resulting resummed expression remains accurate to NLL in the 
variables N and b. 

In this spirit, we make the following approximation, accurate to corrections that are sup- 
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pressed by factors of Qt/pt- 



S 



N+X 



4(p 



r 



\Qt? 



X 



r) 



-N-l 



exp <^ (N + l)Q T -p T /p 



l + O 



( NQ\ 

V PT 



(9) 



Notice that the exponent reaches order unity at just those values of Qt for which the profile func- 
tion begins to decrease. Replacing the singular power dependence on Qt with the exponential, 
we retain the leading behavior at low Qt, but eliminate the kinematic singularity, as desired. 
Again we emphasize that suppression for Qt > pt/N is a reflection of energy conservation. 

Let us now study the effect of the approximation in Eq. (9). Consider for the moment 
N + 1 = — iN with N fixed and real. Using Eq. (7), we may then replace the second line of 
Eq. (4) according to 



C 



ab—t'yc 



A? 

(2tt 



N+X 



c ab ^ ic (4) 



T 

2 \-N-l 



P ab [N,Q T ,^,H 

Xt 



d 2 b exp [E ab ^ c (N, b, Q, /i)] 



d 2 Q 



(27T) 



T e -iAfQ T -p T /p2,-ib-Q T 



(10) 



Here we have extended the Qt integral to infinity. The integral may then be performed, and gives 
5^ (b + A Pt I 'p't) • Using this delta function to perform the b integral in (10), and inserting the 
result back into Eq. (4), we find 



, (resum) 
3 ao AB->-yX 

Pt — — — 

dp T 



E 

ab 



I 

Pt 
8ttS 2 



dN ~ 

^— <Pa/A{N, fl)(/) b/B (N, H) 



N \M ab (x 2 T )\ 2 



x c ah ^ c (4) 



-N-X 



exp 



E, 



ab— >7c 



1 

N, -% ,Q,fi 

Pt 



(11) 



This expression for the direct photon cross section is similar to the result for pure threshold 
resummation, except for the additional b dependence, which has become dependence on the 
combination (N + 1)/pt in the exponent. Although we have derived this form for imaginary 
values of N + 1, it can be analytically continued to any A r , and we use (11) as the result of 
the Qt integral in the high-x-r jointly resummed cross section. Recoil is self-consistently taken 
into account through the exponential in (9), which is accurate up to power corrections as shown. 
There are no kinematic singularities at large Qt- For Qt competitive with p?, of course, the 
approximations we have made fail, but in this region the profile function is small. 



2.4 Matching 

Matching is now straightforward for the stabilized cross section, Eq. (11), and can be handled 
as for the threshold-resummed cross section. We simply expand the exponents to NLO (for 
example) in terms of a s (pT), and replace these approximate expressions with the exact hard 
scattering cross sections at that order. 

We emphasize that within our new treatment we have been able to perform both the Qt 
and the b integrals. This is a great advantage for phenomenological applications, since now 
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the evaluation of the cross section is technically equivalent to that of a standard threshold- 
resummed one. In fact, as we show below, Eq. (11) is identical to normal threshold resummation 
to NLL, but differs at NNLL through a well-defined set of terms that can be identified uniquely 
as recoil effects. The fact that the NLL threshold logarithms are unchanged by recoil is an 
important consistency check of our approach because these logarithms are uniquely specified in 
the perturbative single-inclusive cross section. 

The final resummed cross section thus has a form that is closely related to matched threshold 
resummation [41]: 



7 NLO 7 resl 2 

3 da A B 3 aa A B 



• (12) 



where c^(N) is the moment of the lowest-order cross section, 

PT 1 (AT \1 (AT \ f\,~2 f~2\N \M ab (x^)\ 2 



**W = 8^2 ^ /A {N,^ h/B {N^) / dxl (4) JV ' , (13) 

O^ti JO a/1 — 



and where the final terms in (12) express our matching to the fixed order (NLO, 0(a 2 )) cross 
section ° & AB by taking out the 0(a 2 s ) expansion of the perturbative part of the resummed 

o rj T6s| 2 

cross section, p'j, 



3 Perturbative and Nonperturbative Exponents 



3.1 Resummed perturbative recoil 

To clarify the relationship between joint and threshold resummation and the implications of our 
new treatment of recoil, we review the N and b dependence of the resummed exponent at NLL 
found in [14]. To all orders, the NLL initial-state logarithms in iV and b are generated from an 
integral, derived using the eikonal nature of soft gluon emission, that extends down to zero scale 
in the running coupling. It may be written in a convenient form as 



E l £(N,b,Q,» = Q) 



(14) 



JO K T t : _ >. 



J (bk T ) K { 



2Nk 



( Q 



in 



Nk 7 

~Q 



where Q = 2pT (see Eq. (3)) is the minimal center of mass energy of the partonic subprocess. 
Here and below, we define 



N = Ne yE . 

The anomalous dimensions A a (a s ) have the familiar expansion A a (a s 

= C a 



(15) 



£ n (a s /7r) n with 



= l -C a K=-C a 

a 2 2 



G, 



67 
18 



7l~ 

~6 



10 



T R N f 



(16) 
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where C q = Cf for quarks and C g = Ca for gluons. The presence of the Bessel function 
Ko(2NkT/Q) reflects the conservation of energy that must be imposed to resum threshold and 
kx enhancements simultaneously. For the analogous exponent in kx resummation, the function 
K (2Nk T /Q) in Eq. (14) is replaced by -ln(fc T /Q) and \n(Nk T /Q) by \n(k T /Q), and the k T 
integral produces logarithms of bQ for any b > 1/Q. In joint resummation, however, b must be 
greater than N/Q to produce logarithms. As a result, for N — > oo the profile function in Qt 
space decreases once Qt > Q/N, as discussed in Sec. 2.2 above. 

Starting from Eq. (14) we isolate the effect of perturbative recoil by separating it from the 
corresponding exponent for threshold resummation. Since threshold resummation is already ac- 
curate to NLL in the transform variable N [14], for consistency recoil must appear first at the 
next logarithmic order, and it does. As we shall see, however, its influence on the 1PI cross 
section need not be negligible in perturbation theory. In addition, the integral over the anoma- 
lous dimension v4(a s (/cy)) through the infrared region suggests a specific set of nonperturbative 
corrections, whose effects we will also study. For initial-state radiation, the form of contributions 
beyond NLL accuracy is given in [14]. Threshold logarithms associated with final-state interac- 
tions beyond NLL will be the subject of a separate investigation. We will argue that they respect 
the pattern for power corrections found here. 

As directed by Eq. (11), we now set b = —i(N + 1)/pt in Eq. (14), noting the Bessel function 
relation Jq{iz) = Iq(z). We then reorganize the equation as 



Ef b (N,b 
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Pt 



;Q,fJ- = Q) 



4P T dk\ 
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A (a s (k T )) 
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(N,pr) + 5E ab}rec {N,p T ) 



(N+l)k 7 
Pt 



Nk 7 
~2pr 



- 1 



Kn 



Nk 



T 



Pt 



(17) 



where we have again used Q = 2p?- We have identified the first term on the right-hand-side of 
Eq. (17) with the exponent for threshold resummation for initial-state logarithmic behavior in 
N [14]. 

The second term on the right side of (17) is the recoil correction. It now amounts simply 
to an ^-dependent correction to the threshold-resummed cross section. As required by the 
self-consistency of NLL threshold resummation, this expression is free of NLL logarithms in N, 
because for small arguments z, 



I (z) ~ 1 + 



Kq(z) ~ — In 







-i) 


z 2 


2 









On the other hand, as z = Nkx/pT becomes large with Re(z) > 0, Io(z) increases as e z /\ / 2nz, 
while K (z) decreases as e~ z / (2z/n), so that 



Uz)K n (z) 



— (Re(z) > 0) . 



(19) 



At fixed coupling, and replacing N + 1 by N in Jo in Eq. (17), the net result is a convergent, 
iV-independent integral, equal to (C a + C&) x (a s /2n) ((2), a modest but still significant contri- 
bution in the exponent. 
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For large values of A, we can readily estimate the effect of the running coupling to perturbative 
recoil, by noting that the combination (Jo — l)-Ko becomes sharply peaked near kx = Pt/N, with 
a width that is asymptotically negligible compared to the scale on which the coupling runs. As 
a result, to NNLL, we may isolate perturbative recoil including the running of the coupling by 
the expression, 

r 771NNLL i n | n \ a siM ) Tl^ ) / on \ 
° E ab,rec = i C a + C h ) — Q{2) . (20) 

We will use this expression below to estimate the effects of perturbative recoil. 



3.2 Nonperturbative corrections from threshold 

The recoil exponent 5 E abjTec (N,p T ) in Eq. (17) provides an estimate of perturbative recoil, and 
is also a guide to nonperturbative power corrections. The most basic observation about these 
corrections is that they factorize and exponentiate, in much the same manner as for event shapes 
in e + e~ annihilation and for the transverse momentum distributions of electroweak boson pro- 
duction. This follows from the form of the resummed exponent, in which the entire dependence 
on the running coupling is through a single, integrated scale, kx- We emphasize that a similar 
result holds for the full eikonal exponent to all logarithmic order. Indeed, the same underlying 
nonperturbative parameters that appear in Drell-Yan cross sections will appear in power cor- 
rections to direct photon cross sections. As noted in [41], power corrections from threshold and 
transverse momentum resummations are separately additive in the exponent. As we shall see, 
this leads to an extra power correction compared to estimates based on transverse momentum 
resummations alone [12, 13]. 

Using the additivity of the nonperturbative corrections, we write for the full exponent 



E 



aft— >7C 



(N,p T ) 



E. 
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PT 

abc 



PT 
abc 
PT 

abc, thr 



5E 



np 



i_ r t-iNNLL 



J ab, rcc ) 



(21) 



where 5 E™£ accounts for nonperturbative contributions from low scales in kr, of order Aqcd- The 
full perturbative threshold exponent at NLL, E^ thr , with initial- and final-state contributions, 
was derived in [34, 35, 36, 43]. As noted above, in this study we derive nonperturbative and 
recoil corrections associated with initial-state radiation only. 

For small to moderate values of N, the integral in E^(N,pt) is perturbatively dominated. 
Nonperturbative corrections are generated by treating Nkx/pr as a small parameter in both the 
threshold and recoil exponents of Eq. (17). Expanding the integrands of both 5 E a b, rcc (N,pT) 
and E 1 ^ thr in Eq. (17) for small fcf., we parameterize the resulting l/p% terms as 



^ EH 



(N + l) 2 + A^ 2 
4p 2 T 



E 

i=a,b 



A^ + A^ln 



2pr 



A 2 
4p 2 



\ X ' (Nk/ P t < 1) , (22) 



rri 

1 i=a,b 



where the A 2 terms come from the threshold (Kq) integral in Eq. (17), while the (A + l) 2 
term is from the recoil (Iq) term. The logarithm in both cases arises from the expansion of the 
function K (2NkT/pT), and, as noted above, its presence can be traced to the imposition of 
energy conservation in joint resummation. The constants A™' n in Eq. (22) are interpreted as the 
nonperturbative content of moments of the running coupling [3, 4, 5], with indices in a notation 
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inspired by [4]. More specifically, these are moments of the anomalous dimensions A a {a a (kr)), 
[6] 



\7 n = £ dk 2 T (k 2 T ) m ~ l A a (a s (k T )) ln n (^ N 

= 9± £ dk 2 T (4) m " 1 a s (k T ) In" (^) +.... (23) 



In Eq. (23), the upper limit k for the k\ integral, which also appears as the scale in logarithms of 
Pt in Eq. (22), is a factorization scale. To isolate a truly nonperturbative coupling, as in Ref. [4], 
we could subtract perturbative contributions to the A's to the order corresponding to our level of 
resummation. Since this process does not change the and N dependence of the expressions, 
and because the nonperturbative parameters appear in the same manner here as in electroweak 
annihilation [12], it is not necessary to provide such an analysis for our purposes. 

As in the case of electroweak bosons [10, 11, 44], and in contrast to event shapes [2, 3, 5, 7], 
only even powers result from the expansion of the Bessel functions in (17), the first of which has 
been displayed in Eq. (22). For N not too large, that is, for AAqcd "C Pt we expect only one or 
two power corrections to be significant, but for larger N, the resummed cross section should be 
supplemented by a function with a more general TV- dependence [6, 7]. In the following section, 
we will study the phenomenological implications of such dependence. 



3.3 The full exponent 

Summarizing our results so far, the full exponent is the sum of a perturbative threshold exponent, 
perturbative recoil and nonperturbative corrections, 

E ab ^ c (iV, z^-tl, Q, p\ = E^ CjthT (N, p T ) + 5 E™£(N, p T ) + 5 E^ , (24) 

where the nonperturbative exponent SE^ is given in (22) above and the NNLL recoil correction 
5 -E^ N r cc by (20). E a b t tbi(N, px) is the full exponent for threshold resummation in prompt-photon 
production, including initial and final state contributions [34, 35, 36, 43] (see Appendix A). 



4 Phenomenology of Power Corrections 

The expressions derived above provide useful information on the phenomenology of power cor- 
rections associated with soft gluon emission. First, pr dependence enters through even powers, 
with a leading nonperturbative coefficient that is identical to that encountered as the coefficient 
of b 2 in electroweak annihilation. For comparison, the latter may be written in terms of the same 
parameters A™'™ as 

5 E*? (Drell - Yan) = ~ 

i=q,q 

In contrast, for the single-particle inclusive cross section in joint resummation the nonperturbative 
corrections in Eq. (22) possess highly nontrivial A-dependence, from recoil directly, as well as 



A, 1 ' 1 + A, 1 ' In 



Q 



K 



(25) 
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from threshold resummation. This implies that these power corrections inherit nontrivial energy 
dependence, and we may expect their effects to change with the overall energy. 

To see the qualitative energy- dependence implied by the nonperturbative exponents derived 
above, we note that the Mellin moment N and the variable \nx T are in a conjugate relationship, 
exhibited in the inverse transform (12), 



1 



N "^3- < 26 > 



Identifying these quantities in the nonperturbative power correction of Eq. (22), and recalling 
that xt = 2px/VS, we immediately see that the nonperturbative exponent is suppressed not 
only by a power of px, but also by a power of In S at fixed pt, 



\1,0 I \i,o 

SKI - a tZ^ ^ 



4p T In 



2pT In 



(27) 



Even though this expression is eventually to be convoluted with the hard scattering cross sections 
and the parton distributions, we may conclude that at fixed pt, the importance of power correc- 
tions will decrease as v^S* increases. At the same time, as xx approaches unity, the coefficient 
of 1/pt diverges, and the nominal power correction may dominate at the edge of phase space. 
Even as this coefficient diverges, however, the logarithm in the numerator eventually changes 
sign, so that for xt close enough to unity the enhancement becomes a suppression. This is not 
an accident, because the presence of the factor of N in the logarithm reflects energy conservation, 
which is respected by joint resummation. To give a realistic estimate of the behavior of these 
nonperturbative corrections, we return to moment space. 

For the dominant form of the nonperturbative exponent at moderate N and px, we are guided 
by Eq. (22). For larger N, however, all power corrections in N may become relevant. To account 
for this, we introduce a function of N/px that generalizes this expression. We will refer to this as 
a shape function by analogy to the discussion of [6]. We tailor the ^-dependence to the behavior 
of the Bessel functions of Eq. (17) for large and small values of their arguments, which both 
depend on the combination N/px- Matching to the small- N behavior of the Bessel functions in 
(18) and to the large-iV behavior in (19), we modify Eq. (22), 



C a + C b (N+l) 2 + N 2 + 



2pj_ 

NQ 



7T Ap T A Q N ^ 

\ PT 

Here the scale fio is a parameter of dimension mass squared, which can be thought of as the 
integral of A a (a s {kx)){n /C a ) over k\ with unit weight in (23). The overall factors of color 
charges in Eq. (28) reflect the proportionality of the coefficients A a to C a . The parameter Qq is 
a scale whose value accounts for non-logarithmic terms. If A 1,1 = in (23), we may identify Qo 
with the scale k in Eq. (23). This is the result found in Ref. [4] with n = 2 GeV. 

To estimate the impact of these nonperturbative corrections for large xt, we are aided by our 
experience with electroweak boson production and with perturbative resummation. As can be 
seen from Eq. (23), the nonperturbative parameters in our approach are related to moments of the 
strong coupling, which suggests some form of universality for them. In a study of Z production 
at the Tevatron we estimated nonperturbative effects; the value obtained may be translated into 
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[(Aj-° + Aj'°)/4] ln(m z /«) « 0.8 GeV 2 in Eq. (28), or Aj'° ~ Aj-'° ~ 0.4 GeV 2 , consistent with the 
result quoted in [4] for k = 2 GeV. This implies for the parameter /x in (28), 

7rA 1,0 

/i (Drell - Yan) = — ^ 1 GeV 2 . (29) 

C F 

The nonperturbative exponent (J-E 1 ^ in Eq. (28) with this value of /io is our best estimate of the 
iV-dependent exponent at large N. Substituted into the full exponent (24) it provides a measure 
of power corrections at large xt- 

To illustrate the influence of these power corrections, we compute the ratio of the cross 
section with threshold resummation plus the nonperturbative term (28) to the cross section 
with threshold resummation alone. We do this for several cases that are directly relevant for 
comparison with experiments: for p Be scattering with fixed-target beam energies E = 530 GeV 
and E = 800 GeV (E706 [22]), for pp and pp scattering with beam energy E = 315 GeV 
(UA6 [21]), and for pp scattering at y/s = 63 GeV (R806 [46]). We use the GRV set of parton 
distributions from Ref. [45]. We normalize all our results to the threshold resummed cross section, 
i.e., Eq. (12) with 5 E™^ = 5 E np = 0. This is advantageous because the dependence on the 
factorization and renormalization scales is small in this cross section, and because our recoil and 
nonperturbative corrections have been defined relatively to it. Note that as implied by Eq. (4) 
our resummation is done for the cross section integrated over all photon rapidities. In principle, 
we should account for the finite ranges of rapidity covered in the various experiments, which 
could be done using the techniques developed in Ref. [36]. However, as implied by the results 
of [36], the dependence on rapidity will be very weak in the ratios we consider here and can be 
neglected for simplicity. Our results are always matched to the NLO cross section as described 
after Eq. (12). We do not take into account a photon fragmentation contribution to the cross 
section. 

Shown in Fig. 1 are results for the energies discussed above, as functions of Xt- The en- 
hancements exhibited in the figure are both small and for the most part relatively flat. We have 
included values of xt far from unity, especially for R806, to illustrate the point that these power 
corrections decrease with energy. For RHIC and CDF energies, the effects of (28) are practically 
negligible, of the order of just a few tenths of a percent. We emphasize, however, that this result 
applies only to extrapolations to small xt of the expressions derived for xt — > 1, and is not 
necessarily representative of the true behavior of the cross section at low xt- 

The moderation of the xt dependence of the cross section at large xt in Fig. 1 associated 
with energy conservation is illustrated by comparison to Fig. 2, which shows the analogous ratios 
when the nonperturbative coefficient is allowed to reflect the b 2 (N + l) 2 jp\ dependence that 
is characteristic of kj- resummation, starting from Eq. (25), rather than (28). The shape function 
then has the same overall quadratic iV-dependence as (28), but lacks the iV-dependence in the 
logarithm and the denominator that reflects the influence of the K function in (17). We thus 
have 



& E ab = »o Z ZZ2 ln ( — ) ( Fl §- 2 ) • ( 30 ) 



C a + C b (iV + 1) 2 (2p 7 
ix Ap\ 

Relative to Fig. (1), these curves show both strong enhancements and marked upturns toward 
increasing xt- 

To complete this discussion, we consider two additional variations of the cross sections com- 
puted with Eq. (28). So far, we have ignored the term in (24) associated with perturbative recoil, 
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Figure 1: Ratios of direct photon cross sections computed with threshold resummation and 
nonperturbative shape function (28) to cross sections with threshold resummation only. The 
curves are given as functions of xt for kinematics relevant for comparison to fixed-target and 
ISR experiments (see text). 



Eq. (20). It is probable that the incorporation of recoil at NNLL would affect the values of the 
nonperturbative exponents. Indeed, since both NNLL recoil and power corrections are derived 
from the same starting expression, Eq. (14), there is a serious potential for double counting. On 
the one hand, for small values of iVQo/pr the recoil integral in Eq. (14) is dominated by k? that 
are outside the soft region kx < Qo- On the other hand, once N > Pt/Qo, the integration region 
that gives rise to the result (20) overlaps the power corrections almost entirely. Nevertheless, it 
is interesting to test the influence of the corrections suggested by Eq. (20). To correct for double 
counting, at least partially, we use a modified estimate for recoil, which has the property that 
for small N it approaches (20), while it vanishes for large NQq/pt, 



(C a + C b ) 



a,(4^/iV 2 ) C(2) 



7T 



C(2) V 2pt 



In ( 1 + e 1 / 2 ^ |er 



NQ 



1 + e V 2 -7a 



NQp 
C(2) PT 



(31) 



For N fixed, this is a leading power contribution, with power-suppressed corrections, which, 
however, conspire to cancel the leading term when NQ /pt 1. Figure 3 shows the same sets 
of ratios as in Fig. 1, including now 

SE abf cc in addition to 5E^, Eq. (28). To avoid double 
counting with the C ab ^ lc coefficients in the cross section, Eq. (11), we subtract the leading term 
(C a + Cf,)a s (4p^)£(2)/27r from the latter. We see a substantial increase compared to the pure 
power corrections, in addition to a moderate slope toward large xt- We do not take the level of 
this enhancement too literally, given our rough treatment of double-counting, but conclude that 
it does demonstrate the possible importance of nonleading logarithms and their interplay with 
the magnitudes of the parameters of power corrections. 
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Figure 2: Same as Fig. 1, but for the /^-inspired nonperturbative shape function (30). 

Finally, we illustrate the possible influence of terms that are nonleading by a power in N. 
Although we have derived Eq. (11) only for xt — *■ 1, the form is sufficiently general that it can 
be extrapolated to any value of xt- Clearly, as we leave the kinematic regions where large N 
dominates, terms that are nonleading by powers are expected to become more and more impor- 
tant. Indeed, nonleading terms may be generated from the low-scale limit of partonic evolution. 
Thinking of the pervasive upturn of experimental cross sections relative to NLO noted long ago 
[25], we assume a phenomenological parameterization for the iV-dependent nonperturbative ex- 
ponents that behaves as N jp\ for large N, is related to the splitting functions, and enhances 
the cross section at low xt- The simplest ansatz of this sort is the following modification of the 
quark-gluon exponent (only), 



The parameter Hi is defined by analogy to fi in (28) and (22), but does not have a direct or 
indirect interpretation in terms of resummed perturbation theory. In Fig. 4 we show the same 
ratios, but now computed with the modified shape functions (32), choosing fii = /i = 1 GeV 2 . 
These ratios indeed show a noticeable upturn toward small xt- We observe, however, that the 
magnitudes of the enhancements are nowhere near those necessary to describe the \ow-xt direct 
photon data, especially of E706 [22]. Since we are now considering terms that are subleading at 
large N, we also make exploratory calculations at higher energies, relevant to comparisons with 
the collider experiments at Tevatron (-y/s = 1800 GeV) and RHIC (^/s = 200 GeV). As one can 
see, rather sharp upturns at pr ^ 5 GeV are a distinct possibility here, if our ansatz in Eq. (32) 
is realistic. We finally note, without claims of physics significance, that it is possible to provide a 
qualitatively successful (x 2 per degree of freedom approximately 1.5) "global" fit of direct photon 



5E i9q) 

w np 
w np 




" np 



(32) 
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Figure 3: Same as Fig. 1, but including the subtracted NNLL recoil exponent (31). 

data from E706, UA6, R806, and even CDF, with the ansatz (32), but only for values of fi\ in the 
range of 10 GeV 2 , which implies ^\Ca/^ ~ 2 GeV 2 . The origin of such a large scale, occuring 
as roughly 2GeV 2 /p T , is at the least not obvious. On the other hand, it could be simply an 
artifact of using Eq. (11) outside the region where the exponent E(N,b) in (17), evaluated at 
b = —i(N + 1)/pti has a straightforward interpretation. Even more serious, however, are the 
potential consequences of such corrections for pion production at collider energies. These issues 
can only be clarified by further work. 

5 Conclusions 

We have presented an analysis of recoil and power corrections from initial state radiation in single 
inclusive direct photon cross sections at large xy. In this limit, we resum logarithmic corrections 
in N and simultaneously control logarithmic and power corrections in NQt/pt, where Qt is a 
measure of partonic transverse momentum. Our new treatment avoids any kinematic singularity 
when Qt is large. The resulting expression is equivalent to threshold resummation at NLL in 
perturbation theory, with NNLL recoil effects. We have also shown that we may exponentiate 
power corrections of the form NQt/pt- 

In the large xt region, leading power corrections enter in moment space as powers of (N/pt) 2 , 
with the leading term multiplied by a logarithm of the form \b.(pt/NQ ). We have observed that 
at large xt power corrections are suppressed relative to expectations based on kj- resummation 
alone. This suppression is attributable to phase space restrictions on initial state radiation near 
partonic threshold. This result raises the possibility of a link between a matched fcr resummation 
similar to that of Ref. [32], at relatively low xt, and a joint resummation at large xt- We 
have presented our analysis for initial state radiation, which includes all kj< dependence in joint 
resummation. A detailed discusssion including the role of final state radiation will be given 
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Figure 4: Same as Fig. 1, but for a nonperturbative function with terms nonleading in N as 
given by Eq. (32). We also show results at y^s = 1800 GeV and ^/s = 200 GeV, relevant for 
comparisons with Tevatron and RHIC data. 

elsewhere. 

Looking beyond direct photon production, we anticipate that similar analyses may shed light 
on single hadron and jet production. A simple, but possibly significant observation is that in single 
hadron cross sections, the relevant scale for power corrections associated with partonic threshold 
and transverse momentum is s, the total partonic cm. energy squared. Because s > z~ 2 {\.p\\ 
with z the momentum fraction associated with fragmentation, nonperturbative effects that are 
inverse powers of s are suppressed by factors of z 2 when expressed in terms of p\. Issues such as 
these will be relevant to an effort to tie together perturbative and nonperturbative effects in the 
full range of inclusive hadronic reactions. 
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Appendix 

In this appendix we provide the explicit forms of the exponents E a b^ lc , as given in [34]. Ac- 
cording to Eq. (8) the exponent is split up into pieces associated with initial and final state 
contributions. According to Eq. (17), within our treatment of recoil, the initial-state exponent 



17 



becomes E 1 ^ thl (N,p T ) + E ab ^ cc (N,p T ). One has 



i=a,b 



hf\\) + hY\x^p T ^^ F ) 



_a s (jj? 



(33) 



where 



/4 0) (A) 



-4 



(i) 



2tt6§ 



[2A + (1 - 2A) ln(l - 2A)] 



(34) 



and 



h\ ; (A,2p T ,/i,/i F ) 



4\ 

2ixbl 



^ 2 



ln 2 (l-2A)+2A + ln(l-2A) 



(35) 



+ 



.4 



(2) 



27r6 \ 



A (1) lnfM 
/i 2 



.4 



(i) 



[2A + ln(l - 2A)1 --i-Aln 
nb 



4p 2 T 



For completeness, we have distinguished between the renormalization scale fi and the factorization 
scale /ii?. The A\ are as in Eq. (16), and we have defined 

A = b a s (fi 2 ) In TV , 
11C A -4T R N F 



b 



127T 

17CI - WC A T R N F - QC F T R N F 



247T 2 



For the NNLL exponent E a b, Tec (N,pT) we obtain 



The exponent for the final state reads: 



C(2) 



E™(N,2p T ^) 



a s (/x 2 



/i 0) (A) + /W(A,2pr,A*) + ^(A) 



4i) 



with 



/f(A) = 2fcJ» (A/2) -/if) (A), 
f£\\2pT,y) = 2h£\\/2,2p T ,^2p T )-h£\\,2p T ,p,2p T ) 
A™ In 2 



+ 



41(A) 



7T&o 
7T&0 

cv 



B (l) 

(ln(l - 2A) - ln(l - A)) - ^- ln(l - A) 



nbn 



ln(l - 2A) In 2 , 
ln(l-2A) In 2. 



Here, 



(36) 
(37) 

(38) 

(39) 

(40) 
(41) 
(42) 

(43) 
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Finally, the coefficients C ab ^ c of Eq. (4) read [34]: 




(44) 



(45) 



where 



K 




67 vr 2 
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